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ABSTRACT 

The  research  concentrated  on  one  major  topic  of  investigation:  statistical 

properties  of  inventory  models. 

The  analysis  dealt  with  comparing  the  operating  characteristics  of  several 
Important  stockage  rules  when  the  underlying  demand  distribution  parameters  were 
estimated  from  a limited  sample  of  historical  data  (say,  13,  26,  or  52  observa- 
tions). Particular  attention  was  given  to  the  frequently  observed  case  where 
the  probability  of  zero  demand  occurring  within  any  time  period  was  significantly 
large  (hence,  say,  half  of  the  historical  observations  indicated  zero  demand), 
but  the  positive  levels  of  demand  were  not  necessarily  small.  The  context  of  the 
analyses  was  a large-scale  inventory  system  containing  thousands  of  individual 
items,  and  thus,  the  computational  burden  of  implementing  the  stockage  policies 
was  recognized  in  the  comparison  of  different,  stockage  rules.  The  research 
approach  encompassed  applied  Markovian  probability  analysis  and  computer  simula- 
tion analysis  containing  sophisticated  statistical  autoregressive  time  series 
methods. 

STATISTICAL  PROPERTIES  OF  INVENTORY  MODELS 


Review  of  Work 


Over  the  past  two  decades,  mathematical  analysis  of  inventory  stockage  models 
has  made  great  progress,  and  real-life  implementation  of  inventory  systems  based, 
at  least  in  part,  on  the  results  of  this  modem  research  has  taken  place.  But 
typically  practitioners  have  applied  theoretically  derived  formulas  to  actual 
situations  in  which  the  model's  underlying  assumptions  are  far  from  being  met.  A 
critical  divergence  between  the  model's  assumptions  and  reality  occurs  with  regard 
to  the  amount  of  knowledge  the  systems'  designer  has  about  the  underlying  demand 
distributions.  With  rare  exceptions,  theoretical  inventory  models  assume  that  the 
probability  distribution  of  an  item's  demand  is  completely  specifiable.  In  the 
most  sophisticated  theoretical  inventory  models,  heavy  reliance  is  placed  on  this 
assumption  in  deriving  an  algorithm  for  computing  an  optimal  policy.  Many  other 
theoretical  analyses  make  gross  assumptions  that  permit  stating  inventory  formulas 
which  require  only  a few  parameters  of  the  demand  distribution;  but  even  then, 
these  approximate  analyses  assume  that  the  small  number  of  parameters  are  known 
with  certainty.  A practitioner  utilizing  such  models  typically  employs 
standard  statistical  procedures  to  estimate  the  parameters  in  the  approximation 
formulas,  and  pays  little  if  any  attention  to  modifying  the  estimates  due  to  the 
special  nature  of  the  application.  There  has  been  a significant  knowledge  gap 
concerning  the  nature  of  the  systematic  performance  and  prediction  biases  that  are 
introduced  when  these  statistical  estimates  are  substituted  for  parameters  in  in- 
ventory models. 

The  seriousness  of  the  problem  stems  in  part  from  the  fact  that  in  most  real 
situations,  only  a very  limited  amount  of  past  data  on  demand  are  available,  and 
frequently  such  data  contain  a large  proportion  of  zero-demand  observations.  For 
example,  in  a 26  week  history  of  demand  for  an  item,  perhaps  no  demand  occurred 
during  10  of  the  weeks,  and  in  the  remaining  16  weeks,  the  demand  values  displayed 
wide  variation  (such  as,  30,  100,  7,  245,  etc.).  We  currently  have  no  scientific 
understanding  of  what  stockage  rules  perform  acceptably  well  in  such  an  environment, 
and  how  much  accuracy  we  have  in  our  estimates  of  future  system  performance  based 
on  the  past  data.  The  following  scenario  will  clarify  further  the  context  for  the 
work  supported. 
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Consider  an  Inventory  manager  who  must  design  a system  of  replenishment  rules 
for  the  stockage  of  possibly  thousands  of  items.  Assume  that,  without  too  much 
difficulty  the  manager  can  specify  the  parameters  of  the  criterion  function,  so 
that  the  manager  has  a direct  way  of  determining  whether  one  system  is  better  than 
another,  provided  that  the  manager  can  obtain  the  relevant  associated  operating 
characteristics  for  each  system.  For  example,  assume  that  the  manager  knows  the 
costs  of  holding  inventory,  making  replenishment  decisions,  reviewing  the  stock 
status  of  an  item,  etc.  Also  assume  that  the  manager  is  able  to  articulate  objec- 
tives or  costs  relating  to  stockouts;  to  illustrate,  the  manager  may  specify  a 
cost  attached  to  a unit  of  demand  that  cannot  be  filled,  or  perhaps  the  manager 
may  set  a limit  on  the  probability  of  a stockout,  or  possibly  the  manager  may  tar- 
get fulfilling  a certain  fraction  of  total  demand. 

In  designing  the  system,  the  manager  must  select  a class  of  decision  rules, 
and  preferably,  pick  optimal  policies  from  the  selected  class.  To  illustrate,  the 
complete  economic  preference  function  of  the  manager  ny  imply  that  if  the  manager 
knows  the  probability  distribution  of  demand,  then  an  optimal  policy  is  of  an 
(s,S)  form:  when  inventory  on  hand  plus  on  order  falls  below  s , order  so  that, 
as  a consequence,  inventory  on  hand  plus  on  order  equals  S . The  research  of 
Veinott  and  Wagner  shows  how  to  compute  optimal  values  for  (s,S)  in  this  circum- 
stance, and  the  research  of  Wagner,  O'Hagen,  and  Lundh  tests  out  various  numerical 
approximations  that,  since  the  publication  of  that  paper,  have  turned  out  to  be 
implementable  in  practice.  But  if  the  manager  faces  considerable  uncertainty 
about  each  item's  demand.  dibtfiibution,  because  the  manager  has  only  a limited  num- 
ber of  past  observations  on  demand,  it  is  by  no  means  obvious  as  to  how  the  manager 
should  proceed.  In  fact,  the  decision  process  for  systems  design  becomes  much  more 
complex  than  the  mere  selection  of  parameter  values  for  (s,S). 

First,  the  manager  must  decide  how  much  past  demand  data  to  actually  use.  If 
the  manager  suspects  that  the  underlying  forces  causing  the  demand  for  an  item 
shift  from  time  to  time  the  manager  may  want  to  ignore  data  that  are  older  than 
some  number  of  periods,  which  the  manager  then  must  specify.  Since  the  system  is 
to  operate  over  an  indefinite  future,  the  manager  also  must  determine  how  often  to 
update  the  rules,  that  is,  how  often  to  discard  old  historical  data  and  recompute 
the  replenishment  policies'  parameters.  For  example,  the  manager  may  choose  to 
revise  the  policies  every  six  months  and,  at  each  revision,  utilize  the  past  twelve 
months'  data. 

Second,  since  the  manager  does  not  know  the  form  of  the  demand  distributions, 
the  manager  must  either  guess  at  the  form,  then  statistically  estimate  the  distri- 
bution's parameters  for  each  item  (possibly  using  a prior  distribution),  and  finally 
compute  policies  assuming  that  the  estimated  distributions  are  in  fact  the  true  dis- 
tributions; or  the  manager  may  choose  to  use  an  "approximately  optimal"  form  for  the 
replenishment  policy  (such  as  let  S equal  a uumber  of  mean  demands  plus  a multiple 
of  the  standard  deviation  of  demand,  and  s»S-D  , where  D is  the  familiar  Wilson 
square  root  lot  size  formula  that  requires  knowing  only  the  mean  demand),  and  then 
select  statistical  estimators  for  the  parameters  required  by  this  policy  form.  This 
enumeration  does  not  make  apparent  the  full  design  decision  possibilities  because, 
within  each  option,  the  manager  has  many  alternatives  from  which  to  choose.  To  com- 
plete the  scenario,  suppose  that  the  manager  does  make  a particular  selection  from 
among  all  these  many  options. 
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Third,  Che  manager  must  decide  fae|$OA£  implementing  the  system  vhecher  the 
design  parameters  are  well  set.  For  example,  if  the  manager  lets  the  stockage 
objective  level  S ■ a* (mean  demand)  + ba  (standard  deviation  of  demand),  the 
manager  must  decide  on  the  appropriate  values  for  a and  b . Since  the  manager 
realizes  that  there  may  have  been  a compounding  approximation  on  approximation, 
the  manager  usually  hesitates  to  rely  solely  on  theoretical  probability  distri- 
bution analysis  that  is  based  on  complete  knowledge  of  the  demand  distribution. 

Most  likely,  the  manager  will  employ  retrospective  simulation  (that  is,  will 
reuse  the  limited  amount  of  past  data  from  which  the  manager  estimated  the  mean 
demand  and  the  standard  deviation  of  demand)  to  adjust  the  values  of  a and  b 
so  that  the  resultant  6taAAAtLcaJL  eAtimaA.il  of  the  criterion  function  is  at  an 
optimal  value.  (Actually,  the  manager  probabily  will  set  the  same  values  for  a 
and  b over  a large  number  of  items,  and  so  will  perform  the  adjustment  process 
via  retrospective  simulation  by  using  an  aggregate  objective  function.)  The  man- 
ager also  will  estimate  several  operating  characteristics  of  importance,  such  as 
average  inventory  on  hand,  the  probability  of  a stockout,  the  fraction  of  demand 
filled,  etc. 

To  summarize  the  scenario,  we  see  that  the  system's  designer  must  select  in 
concert  the  number  of  historical  observations  to  use,  the  frequencey  for  repeating 
the  reestimation  process,  the  form  of  the  replenishment  rule,  the  statistical 
estimators  to  produce  the  demand  parameters  required  by  the  rule,  and  the  design 
parameters  of  the  rule.  The  manager  makes  all  of  these  choices  based,  at  least  in 
part,  on  simulating  how  the  proposed  system  would  have  performed  in  the  past  (and 
in  doing  so,  the  manager  typically  uses  the  same  limited  data  for  both  estimation 
and  prediction).  In  all  the  discussion  to  follow,  we  employ  the  term  HJ&Aojn' & dihAgn 
to  mean  the  entire  composite  of  these  many  choices. 

The  nature  of  the  supported  research  was  to  examine: 

1.  How  good  are  the  manager's  statistical  estimates  of  the  system's  future 
performance  for  selected  choices  of  the  system's  design?  Similarly,  how  good  are 
theoretically  based  estimates  of  future  performance? 

2.  Under  what  circumstances  do  certain  "approximately  optimal"  policies 
perform  better  than  others?  Are  some  policies  statistically  more  robust  than  others? 

3.  How  do  the  statistical  estimates  of  the  system's  future  performance,  as 
well  as  the  performance  characteristics  themselves,  depend  on  choice  parameters  such 
as  the  number  of  historical  observations  to  use  and  the  frequency  of  reestimating 
the  policy's  demand  parameters? 

Special  attention  was  given  to  systems  in  which  the  demand  distributions  dis- 
played a high  frequency  of  zero  demand  (the  positive  levels  of  demand,  however,  did 
not  need  to  be  anywhere  near  zero) . 

Except  under  extremely  special  assumptions,  it  was  not  possiple  to  derive  com- 
putationally practical  formulas  for  the  exact  statistical  distributions  associated 
with  the  above  questions.  The  best  hope  for  increasing  our  ability  to  answer  the 
research  questions  above  in  a wide  variety  of  situations  was  to  derive  more  workable 
numerical  approximations.  With  that  objective  in  mind,  we  built:  under  the  re- 
search grants  a computer  simulation  model  that  yielded  the  actual  operating  charac- 
teristics and  the  corresponding  properties  of  statistical  estimators  of  these  opera- 
ting characteristics  for  inventory  systems  that  contained  uncertainty  about  the 
demand  distributions.  The  details  of  the  simulation  design  are  discussed  below. 


Exptotutioxy  Simulation  Expe/umenti.  The  selection  of  the  underlying  probability 
distributions  and  other  parameters  used  in  the  research  was  influenced  by  the  earlier 
study  of  Wagner,  O'Hagen,  and  Lundh.  The  demand  distributions  were  both  Poisson  and 
Negative  Binomial  (variance/mean  - 3,  9),  as  well  as  distributions  that  were  prob- 
ability mixtures  of  zero  demand  and  a standard  Negative  Binomial  distribution 
(variance/mean  - 9).  More  than  144  sets  of  the  economic  and  lead-time  parameters 
were  investigated.  A summary  of  the  case  selections  is  shown  in  Table  1;  unit 
holding  cost  was  normalized  to  equal  1.  A few  additional  cases  with  higher  mean 
demands  also  were  run. 

Table  1 


Demand 

Distribution 

Mean 

Demand 

Unit  Penalty 
Cost 

Setup 

Cost 

Lead 

Time 

Poisson 

2,4,8,16 

4,9,99 

8,32,64 

0,2,4 

Negative 

Binomial 

(V/p  - 3,  9) 

2,4,8,16 

4,9,99 

8,32,64 

0,2,4 

Research  concentrated  on  investigating  optimal  and  approximately  optimal  (s,S) 
policies  assuming  that  demand  for  each  period  is  independently  and  identically 
distributed.  For  statistically  derived  policies,  the  formulas  in  Wagner,  O'Hagen, 
and  Lundh  were  used,  where  every  J periods  the  demand  history  over  the  past  K 
periods  was  employed  to  update  the  estimates  of  the  mean  and  variance  of  the  demand 
distribution.  The  values  for  (J,K)  that  have  been  tested  are  (13,13),  (13,26), 
and  (26,26).  Each  time  the  policies  were  revised,  forecasts  of  important  operating 
characteristics  were  made  using  restrospective  simulation  over  the  past  L periods. 
The  values  for  (J,  K,  L)  that  were  tested  are  (13,13,13),  (13,26,13),  (13,26,26), 
and  (26,26,26). 

Summary  of  Results 

An  illustration  of  the  simulation  model  output  is  shown  in  Figure  1.  The  case 
was  a Negative  Binomial  distribution  for  demand  with  a mean  of  4 per  period,  a lead 
time  of  4 periods,  a unit  penalty  cost  of  99,  a setup  cost  of  64,  and  a unit  holding 
cost  of  1.  An  optimal  policy  was  (33,57)  under  these  assumptions,  with  expected  cost 
per  period  equal  to  40.50.  [This  policy  was  obtained  in  the  Wagner,  O'llagen,  and 
Lundh  study.]  The  simulations  in  Figs,  la  and  lb  were  run  for  200  cycles,  each  con- 
taining 13  periods;  the  simulations  in  Fig.  Id  were  run  for  200  cycles,  each  contain- 
ing 26  periods.  Simulation  estimates,  and  their  standard  errors,  are  shown  for 
period-end  inventory  level,  stockout  quantity,  stockout  frequency,  replenishment 
quantity,  replenishment  frequency,  and  total  cost  per  period.  The  Normal  approxima- 
tion policy  (32,55)  also  was  simulated — in  this  particular  illustration,  it  outper- 
formed the  optimal  policy,  because  of  the  particular  realized  history  of  random 
demands.  The  bottom  half  of  Figs,  la,  lb,  and  Id  show  the  behavior  of  the  statisti- 
cally based  policies  and  the  corresponding  forecasts.  In  the  (13,13)'  simulation, 
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the  average  cost  per  period  for  the  statistical  policy  vas  51.52  (standard  error 
3.549) , and  the  forecasted  average  cost  per  period  using  13  periods  for  the  retro- 
spective simulations  vas  34.66  (standard  error  .618). 

An  illustrative  tabulation  of  resulting  total  cost  per  period  under  optimal 
(s,S)  policies,  approximately  optimal  policies,  and  statistical  policies,  is  shown 
in  Fig.  2 for  all  cases  of  Negative  Binomial  demand  with  lead  time  4 . Similar 
tabulations  are  available  for  all  the  Poisson  cases  and  the  Negative  Binomial 
cases  with  other  lead  times.  Similar  tabulations  also  are  available  for  other 
operating  characteristics,  including  standard  deviation  of  total  cost  per  period, 
period-end  inventory  level,  backlog  as  a proportion  of  mean  demand,  frequency  of 
backlog,  and  replenishment  frequency. 

An  illustrative  tabulation  of  the  forecast  biases  is  given  in  Fig.  3 for  all 
the  cases  of  Negative  Ninomial  demand  with  lead  time  4 . Average  values  for  the 
statistically  determined  (s,S)  policy  parameters  were  tabulated,  as  illustrated  in 
Fig.  4.  Th®«*e  averages  were  remarkably  close  to  the  optimal  values  for  (s,S). 

Because  of  the  considerable  magnitude  in  many  of  the  standard  errors  of  the 
estimated  mean  values  for  the  operating  characteristics,  the  above  tabulations 
provided  "at  a glance"  ouly  orders  of  magnitude  and  plausible  relationships,  but 
still  left  open  to  question  the  regularity  of  certain  observed  behavior.  Hence, 
another  means  of  summarization  was  used  to  clarify  the  results. 

The  cases  shown  in  Table  1 above  were  aggregated  into  multi-item  inventory 
systems,  namely,  the  72  Poisson  demand  cases  and  the  two  sets  of  72  Negative  Binomial 
cases.  A large  number  of  such  summarizations  for  the  various  operating  character- 
istics were  tabulated,  and  a few  will  be  illustrated  below  to  show  their  format 
and  suggest  their  implications. 

Figures  5a  and  5c,  for  Poisson  and  Negative  Binomial  (variance/mean  - 3)  demand, 
respectively,  tabulated  the  proportion  of  total  cost  associated  with  the  various 
parameter  values  in  each  72  item  system  that  employed  statistical  policies  which 
were  recomputed  each  26  periods  using  26  periods  data.  For  example,  in  the  Poisson 
system,  the  total  system's  cost  per  period  was  2037.8,  of  which  39.9%  was  due  to 
items  with  penalty  cost  99,  56.6%  to  items  with  setup  cost  64,  38.1%  to  items  with 
mean  demand  16,  15%  to  items  with  penalty  cost  99  and  mean  demand  16,  etc.  Figures 
5b  and  5d,  for  Poisson  and  Negative  Binomial  demand,  respectively,  tabulate  the  per- 
centages by  which  the  statistical  system's  total  cost  per  period  exceeded  that  of 
the  optimal  policy.  For  example,  in  the  Negative  Binomial  system,  total  system 
cost  for  statistical  policies  exceeded  that  of  optimal  policies  by  only  5.8%. 

Similar  multi-item  systems  summaries  were  tabulated  for  the  Important  operating 
characteristics  of  interest,  including  period-end  inventory,  backlog  cost,  backlog 
frequency,  replenishment  cost,  and  replenishment  frequency. 

The  forecasts  were  compared  in  a similar  manner.  Sources  of  the  underestimation 
bias  for  a statistical  system  that  was  revised  every  26  weeks  using  26  weeks  of  past 
data  to  set  the  stockage  parameters  and  retrospectively  simulate  is  shown  in  Fig.  6 
for  the  Negative  Binomial  items.  In  this  illustrative  example,  total  cost  per 
period  was  underestimated  by  9.3%  for  the  total  system  comprised  of  72  items.  Simi- 
lar summaries  were  tabulated  for  other  settings  of  the  review  frequency  and  for 
other  operating  characteristics. 
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The  following  reports  summarize  the  results  and  implications  of  all  the  des- 
cribed simulation  experiments: 

STATISTICAL  PROBLEMS  IN  INVENTORY  CONTROL,  Alastair  MacCormick,  December 
1974,  244  pp. 

MULTI-ITEM  INVENTORY  SYSTEM  POLICIES  USING  STATISTICAL  ESTIMATES:  NEGATIVE 

BINOMIAL  DEMANDS,  Arthur  S.  Estey  and  Ronald  L.  Kaufman,  September  1975, 
85  pp. 

VARIANCE  REDUCTION  TECHNIQUES  FOR  AN  INVENTORY  SIMULATION,  Richard  Ehrhardt, 
September  1975,  24  pp. 

COMPUTER  PROGRAMS  FOR  (s,S)  POLICIES  UNDER  INDEPENDENT  OR  FILTERED  DEMANDS, 
Ronald  L.  Kaufman,  (in  preparation),  71  pp. 

MULTI-ITEM  INVENTORY  SYSTEM  POLICIES  USING  STATISTICAL  ESTIMATES:  SPORADIC 

DEMANDS,  Ronald  L.  Kaufman  and  John  G.  Klincevicz,  (in  preparation), 

83  pp. 

THE  POWER  APPROXIMATION,  A NEW  ALGORITHM  FOR  CONTROLLING  INVENTORY  SYSTEMS 
WITH  INCOMPLETE  DEMAND  INFORMATION,  Richard  Ehrhardt,  (in  preparation), 

102  pp. 
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Figure  5a 


Figure  5b 
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Figure  5c 
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